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Motivated by recent experimental results reporting giant coercive fields in Co(II)- 
based molecular magnets we present a theory of hysteresis phenomena based on the 
Glauber stochastic dynamics. Unusual form of hysteresis loops is similar to those 
of found in Co-based quasi-one-dimensional ferrimagnet CoPhOMe at low tempera- 
tures. Temperature dependence of the coercive field has a characteristic form with an 
inflection that may serve as an indicator of the Glauber dynamics in real compounds. 
A relevance of the model for other Co-based molecular magnets is discussed. 

PACS numbers: Valid PACS appear here 

I. INTRODUCTION 

One of the remarkable phenomena found in molecular magnetic materials is the magnetic 
hysteresis similar to those observed in hard magnets.-* 2 - In particularly, magnetic hysteresis 
can be observed in the absence of a long-range magnetic order in zero-dimensional single- 
molecule magnets (SMM) 3 and one- dimensional (ID) single chain magnets (SCM). 4,5 The 
hysteresis phenomena in SMM are affected by large easy-axis magnetic anisotropy and by 
the weak intermolecular interactions.- Although a large easy-axis magnetic anisotropy can 
also be of importance in quasi-lD magnets, an origin of magnetic hysteresis in these systems 
is not yet well understood.- 

Discovery of slow relaxation of magnetization in quasi- ID compound 
Co(hfac) 2 [NIT(C 6 H 4 p-OMe)] (or CoPhOMe)^ evoked interest to the Glauber kinetic 
model suggested for relaxation dynamics in ID Ising ferromagnet.-"^ The source of a 
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strong easy-axis anisotropy are Co (II) ions that enable to treat collective properties of this 
compound by means of ID Ising model. Recently, pronounced hysteresis loops with giant 
coercivity were observed in another Co(II)-based materials, namely [Co(hfac)2NIT-C6H4- 
O-R] 11 , where R=(CH 2 ) 3 CH 3 , and [Co(hfac) 2 ]-BNO* (BNO* is chiral triplet bis(nitroxide) 
radical).— Despite the Glauber dynamics seems to be plausible, one should take into 
account that chains are packed in a 3D structure, and with a cooling ID units undergo a 
phase transition into 3D magnetic array. It means that another explanation of hysteresis 
phenomena related with domain wall (DW) dynamics comes into the play. High coercive 
fields reported in Refs.— ^ at low temperatures, i.e. broad hysteresis loops, support the 
reasonings. Note that considerations of slow magnetic relaxation below 3D ordering in 
the framework of the DW approach turns out to be fruitful for Mn(II)-based quasi-lD 
ferrimagnet s .— ^ 

Being guided this motivation we present a theory of magnetic hysteresis for Ising ferrimag- 
netic chain calculated within the Glauber dynamics approach. An extension of the Glauber 
model for higher spins (more than 1/2) has been done in a modeling of photoinduced re- 
versible magnetization^ 1 ^. Despite equilibrium properties of mixed spin Ising systems have 
attracted a much of attention, a study of nonequilibrium aspects of the model started only 
recently— In the paper we perform the mean field treatment of the Glauber-type stochastic 
dynamics of the mixed (3/2,1) ferrimagnetic chain by considering fluctuations of local fields 
in the spirit of the generalized mean-field theory— Our goal is to verify a validity of the 
model for real Co-based molecular ferrimagnetic chains. The choice of the spin quantum 
numbers is stipulated by a possible prototype of the model, the metal-organic compound 
[Co(hfac) 2 ]-BNO* studied recently by two of us, which embodies Co(II) ions with spin 3/2 
and the chiral triplet biradical ligands BNO* with spin 1. However, the results we obtain 
may be of interest for another Co-based quasi-lD magnets. 

The paper is organized as follows. In Sec. II we obtain dynamic equations of the Ising 
ferrimagnetic {S,a) chain model derived through the Glauber stochastic dynamics in a pres- 
ence of a time-varying magnetic field. In Sec. Ill we present results of hysteresis behavior 
obtained via numerical calculations of these equations for the case 5 = 3/2, a = 1. A depen- 
dence of shape of the hysteresis cycles on the ratio between the magnetic field frequency and 
the spin flip frequency is studied. In Sec. IV we confirm by a Monte Carlo (MC) method 
that a hysteresis loops arise in a static magnetization process. Discussions are relegated to 
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the Conclusion part. 



II. MODEL 



The Hamiltonian of the Ising (3/2,1) ferrimagnetic chain is given by 



(1) 



i \ i i 



where the first term sums interactions between the nearest neighbors with the spin variables 
Oi = 0, ±1 and Si = ±1/2, ±3/2, and J > favors an antiferromagnetic coupling of the 
adjacent sites. The Zeeman term describes interaction of the spins with an oscillating 
magnetic field of the sinusoidal form 



with a frequency lu. The system is in contact with an isothermal heat bath at given tem- 
perature T. 

The relaxation of the interacting Ising system cannot be obtained from its Hamiltonian 
because it eliminates intrinsic spin dynamics (precession in the local field). Nevertheless, 
relaxation phenomena can be described by means of a phenomenological equation which 
specifies the transition rate from one spin configuration to another. In general, we assume 
that transition from one configuration to another involves changing a single spin, i.e. ac- 
cording to Glauber dynamics, -^>^ at a rate of 1/r transitions per unit time. 

According to formalism of Ref.— the master equation for one sublattice can be written 
in suggestion that spins on another sublattice momentarily frozen, i.e. 

d . 

— -P{s}(o"i, cr 2 ,...,(Tj,..., a N ) 



H{t) = H Q cos(ujt) 




3 T'j^Tj 
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dt 
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7^ ^'(Sj ~ * Sj)P{cr}{Sl, S2, . . . , Sj, . . . , S 



"EE ^ ^ P M(Si> & S S N \ (3) 

where P{s}(&i, &2, ■ ■ ■ &j, • • • , &n) is the probability that the system has the spin configuration 
{o"i, cr 2 , . . . <Jj, . . . , cr/v} in the first sublattice leaving the spins {5*} of the second sublattice 
fixed, Uj(aj — > <Jj) is the probability per unit time that the j-th spin changes from the value 
a'j to <jj. The similar notations hold for another sublattice. 

The transition probabilities u)j(a'j — > Oj) and ujj(S'j — > Sj) are imposed by the principle 
of detailed balance. Indeed, the principle requires the probabilities of states P{$} and P{ a y 
to be stationary at equilibrium 

iOj((jj -> a'j)P{s}{ox, o 2 , • • • ,crj, • • • ,0jv) = -> ooO-PjsjOi, 02, . . . ,erj, . . . ,crjv), 

~^ Sj)P{ a y(Si, S*2, . . . , Sj, . . . , Sn) = ojj{Sj — > 5' jl -)P{ -}(S'i, S2, . . . , Sj, . . . , Sn)- 
Provided the probabilities P{s} and P{ CT } are Boltzmann distributions, i.e. 
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where = /? [F - J(5 i _ 1 + S,)] , and ^ = /? [P - J(oj + o i+1 )] , (3 = 1/T. 
In Glauber dynamics the relationships ([6]) are resolved as follows 
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where Q = 1/r is a number of spin changes per unit time. Strictly speaking, Eqs.([7j) are 
relevant for a equilibrium process with a constant magnetic field. Nevertheless, it is reliable 
when the field sweep frequency is much less then that of spin transitions u <C Q. 
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Expectation value of the cr-sublattice magnetization at the moment t is given by 



(8) 



where 



= Yl Yl ' P {S}(^1, ...,<Ti,...,(T N ) 

{S} M 



(9) 



is the probability to find the <7, value. The sum with the prime runs over all a variables 
except that of i-th site. Similar definitions are hold for another sublattice. 
The probability satisfies the following dynamic equation 



dt 



Eexpfo'y, 



!„. 1 / ' 
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which is obtained from Eq.([9]) with an account of Eqs. (J2H3]) and (j7j). The average in the 
right-hand side is determined as follows 



EE 



exploit 
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The analogous equation is derived for another sublattice 



dt 



si 



The system of Eqs. fHOII12p can be treated in the mean-field approximation (MFA) 
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that yields 
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where the Brillouine function B s (x) is introduced 
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FIG. 1: Hysteresis loops of magnetization per cell M = (a) + (S) at different temperatures 9. 
Inset: the coercive field h c vs temperature. 

By taking the uniform arrangement (crj) = (a), (Si) = (S) and, as a consequence, (yj) = 
(3[H{t) — 2J(S)}, (£,•) = /3[H(t) — 2 J (a)] one obtain eventually the dynamical equations 



dt 



d(S) 
dt 



-n 



-n 



(a}-aB a (a(3{H(t)-2J(S)}) 



(15) 



(S)-SB s (Sp{H(t)-2J(<r)}) 

At equilibrium (d/dt — ■> 0) one recover the usual MFA equations from the system. 

In one dimension the mean-field approximation is poor since the local fields yi and £j 
fluctuate strongly from one site to another. To overcome partly the drawback we use a 
generalization of the MFA approach like those used in Ref.— 

Take the following approximation for the averages ( TTTT) 
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i.e. the averages depend on the spin variables at the adjacent sites of another sublattice. 

After substitution Eqs. (ll6HTTj) into Eqs.( iTU|T2"l) one obtain the dynamic equations of the 
generalized mean-field approximation written through the probabilities 

d 
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By determining a time evolution of these quantities one recover a time dependence of ex- 
pectation values for the sublattice magnetizations according to Eq.([H]). 

III. NUMERICAL RESULTS 



Let us now discuss hysteresis phenomenon for (3/2, 1) ferrimagnetic Ising chain obtained 
on the base of Eqs. (1151) and (fl8l) . For numerical calculations it is convenient to rewrite the 
MFA equations in the form 

d(a) 



dr 



d(S) 
df 



■(a)+B l ([h-2(S)]/9) 



■(S) + (3/2)53 (3[/i-2(<7>]/20). 



where f = fit, and h = H/ J, 9 = T / J are reduced field and temperature, respectively. 

Regarding the system fflBl there is a way to bring down a number of differential equa- 
tions from 7 till 5 by using the normalizing condition for probabilities p and definition of 
observables. Indeed, hold four differential equations for the probabilities p(±l), p(±3/2) 



d r-i-i\ i V exp[±(h-Si-i-Si)/0\ , , , , , 

- P (±i) = - P (±i) + £ ^^h-s^-sm ^ l - l)p{ t] ' (19) 



|p(±3/2) 



-M±3/2) + £ expl ^ (ft -^£f 9)1 ^M<W. (20) 
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and add the corresponding equation for the observable (S) 
d(S) ^ 3 



df 

obtained from Eq. fllSI) . The system is supplemented by the pure algebraic relations 

p(0) = l-p(-l)-p(l), 

p(-l/2) = p(3/2) - 2p(-3/2) - (S) + 1/2, 
p(l/2) = (5) - 2p(3/2) +K-3/2) + 1/2 
that yield the leaving probabilities. 



The sums in the right-hand sides of Eqs. (I19II20I) are explicitly given by 

_ exp[o-j(/i - - Si)/ 6} wo ^ _ exp[(7i(/i + 3)/6>]^ o Mi2 



^ 2^exp[cr;(/i- S<_i - Si)/0] Zi[(/i + 3)/0] 



+ ^ Z P 1 1 ^' [?(-3/2)p(3/2) + p(-l/2)p(l/2)] + eX ^f/ 1) 1 / y 1 [p(l/2)V2p(-l/2)p(3/2)] 

for the spin-1 sublattice and 

V- e^[S l {h~a j -a j+l )/e} _ exp^fe + 2)/fl] 2 



for the spin-3/2 one. Here, 

Zi(x) = ex P [ a 'i x \ = 2 cosh(x) + l, Z:i(x) = exp[^x] = 2 cosh(x/2)+2 cosh(3x/2). 

The results of such calculations performed within the Runge-Kutta method are presented 
in Figs. 1-3 for the field frequency u/Q = 10~ 4 (Q = 1). Initial values were chosen to 
correspond either to total saturation or disorder in both sublattices. A stationary regime is 
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approximately reached after one full sweep. Figfj] shows the evolution of hysteresis loops at 
various temperatures. It can be clearly seen that the area of the hysteresis loop monotonically 
decreases with increasing temperature. The hysteresis curves are characterized by well- 
defined steps, and for higher temperature (6 > 0.7) have the form with a loop in the middle 
close to that of observed experimentally in Co(II)-based quasi-lD ferrimagnetic compound 
CoPhOMe at temperatures around > 4.5 K.- The process of thermally activated spin flip 
formation reduces the effective intrinsic coercive field (see inset in FigJT]). Note that the 
temperature dependence of the coercivity exhibits a behavior (a curve with an inflection) 
distinguished from that of predicted by a thermal activation theory by Egamt^ used for 
hard magnetic materials.— When T tends to zero the value of the coercive field becomes 
equal to 2 J, i.e. it is determined by the exchange coupling (not by anisotropy) as expected 
for Ising systems. 

The both analytical approaches give qualitatively same results (see Fig. 2). Note only 
that a more detailed account of fluctuations within the generalized MFA squeezes the area 
of the hysteresis loops. 

To elucidate physics standing for this hysteresis behavior, magnetizations of the both 
sublattices are picked out in Fig. 3. The results resemble a spin-orientational phase transition 
in the Heisenberg antiferromagnet and are interpreted as follows. At strong fields both 
sublattices are polarized along the field. With its decreasing the bigger spins retain their 
directions while smaller spins change smoothly their alignment into opposite to gain in an 
exchange energy. During further demagnetization process, when h < 0, the bigger spins 
reorient to be again arranged along the field. This causes a sharp spin reorientation in 
another sublattice due to the exchange coupling which is stronger than the corresponding 
Zeeman interaction. The rapid change is accompanied by the side effect which stands out 
more noticeably in the MFA calculations, namely, an appearance of "whiskers" in the entire 
hysteresis curve at low temperatures. The process is completed by a gradual saturation 
of small-spin sublattice magnetization along the applied field. The back sweep goes on 
similarly. 

We also found out how the magnetization process depends on the realtionship between the 
frequency Q (a number of spin transitions per unit time) and the magnetic field frequency 
to. We plotted the corresponding curves in Fig. 4 where two cases are presented, namely, 
the quasistatic regime uo/fl = 10 -4 and the regime when the discrepancy between both the 
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FIG. 2: Comparison of hysteresis curves calculated within MFA (solid) and generalized MFA 
(dashed) approaches. 

frequencies are not so drastic uj/Q — 10 _1 . As well we show in the Figure the pure static 
magnetization process which almost coincides with that of the quasistatic regime. One see 
that in the case uj/VL = 10 _1 the hysteresis loop transforms into the narrow S-like form 
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FIG. 3: Hysteresis loops for the sublattice (a) (big dotted), (S) (small dotted) and the entire M 
(solid) magnetizations (Q = 0.4). 

similar to those found experimentally in the compound CoPhOMe at 2.0 K and 3.5 K. 
Physically, the increasing the ratio oj/Vt at a fixed magnetic field frequency means that a 
spin dynamics governed by Q slows down. 

IV. CONFIRMATION BY MONTE CARLO SIMULATIONS 

In the previous section we have shown that in the quasistatic regime u <C Q the both 
analytical models demonstrate an appearance of hysteresis loops of a peculiar form in the 
middle of magnetization curves. However, we note that spin fluctuations taken into account 
in the framework the generalized MFA are strongly restricted since they hold a translational 
symmetry of the chain, i.e. they are the same within an each elemenatry cell. In order to 
check whether the hysteresis phenomena are stable against fluctuations in common case, we 
study the static magnetization process by a Monte Carlo method. 

We apply standard importance sampling methods^ 1 ^ to simulate the Hamiltonian given 
by Eq. ([1]). Periodic boundary conditions on iV = 64, 256 chains were imposed and configu- 
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FIG. 4: Evolution of magnetization curves for different frequency regimes (6 = 0.2): (a) MFA; (b) 
generalized MFA 
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rations were generated by sequentially traversing the sublattices and making single-spin flip 
attempts. The flips are accepted or rejected according to a heat-bath algorithm. Our data 
were generated with 10 4 Monte Carlo steps per spin in the chain after 10 3 warming steps 
per spin. We checked that the size effects do not substantially effect the results (Fig. 5). 
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FIG. 5: MC magnetization curves for different chain lengths. 

Qualitatively, the same features of the hysteresis loops as obtained in MFA and the 
generalized MFA are also found in MC. Namely, the hysteresis loop of an almost ideal 
rectangle form appears in the middle of the magnetization curve at low temperatues. The 
temperature evolution of the hysteresis curves reproduces qualiatatively that of predicted 
by the analytical treatments, however, the MC hysteresis loop narrows more rapidly with 
an increasing of temperature (Fig. 6). 

V. CONCLUSIONS 

We analyse a magnetic hysteresis of Co-based quasi- ID ferrimagnetic magnets within 
the model of the mixed spin Ising chain. By using a Glauber dynamics approach we build 
hysteresis loops that come up when a sinusoidal magnetic field is applied. We found that the 
unusual shape of the calculated hysteresis cycles coincide with those found experimentally 
in the CoPhOMe ferrimagnet at low temperatures. However, another Co-based molecular 
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chainsii^ demonstrate a hysteresis behavior that is scarcely agreed with the analytical 
treatment, i.e. these materials behave similar to very hard magnets with a high coercive 
field and broad hysteresis loops. 

It is likely that a reason behind the hardness of the materials is interchain interactions 
that result in 3D ordering at low temperatures, and, as consequence, change a character 
of elementary excitations. Namely, random spin flips throughout the Ising chain, that is a 
feature of the Glauber dynamics, are substituted for spin flips in the vicinity of domain walls 
separating regions of opposite magnetizations that causes their displacements. Then the 
demagnetization process is related with the thermally activated DW motion. Similar ideas 
have been recently argued in Ref.-, where the Mn(III)-based systems^ 3 - 1 ^ were suggested to 
find out a role of easy axis anisotropy in hysteresis phenomena. 
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FIG. 6: Temperature evolution of MC hysteresis loops in a comparison with the mean field (MF) 
and the generalized mean field (GMF) calculations (N = 256). 



